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DESCRIPTION 

METHOD AND APPARATUS FOR PREDICTING BENDING LIFE SP ANS OF 
ELECTRIC WIRES AND/OR WIRE PROTECTING MEMBERS^_ INgU^^ 
VIB RATIONS, AND RECORDING MEDIUM STORING PRO GRAM 

5 

Technical Field 

The present invention relates to a method and apparatus 
for predicting bending life spans of a plurality of electric 
wires clamped at least two points to a predetermined part and/or 
1 0 their wire protecting members , the bending induced by vibrations , 
and a program for predicting the same. 



Background Art 

Usually, in a vehicle or the like, a plurality of electric 

15 devices are electrically connected to one another by a plurality 
of electric wires (referred to frequently and simply as wires) . 
The wires are bundled, by an insulation locking band, tape or 
the like, into a single bundle as so-called wire harness which 
is extendedly arranged at predetermined locations of the 

20 automobile or the like. In some cases, the wires are protected 
by a wire protecting member such as a grommet are arranged at 
a predetermined part. In other cases, the wires are arranged 
at a predetermined part of the vehicle or the like, while not 
bundled into one line and/or not protectedby the wire protecting 

25 members. 
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The vehicle is always placed in such an environment that 
it is always vibrated due to engine drive and the like. In 
particular, the electric wires distributed inside the engine 
room directly receiving the engine vibrations, are repeatedly 
bent and deformed, and finally disconnect due to the vibrations . 
Experimental results of such are reported- In this respect, 
it is of particular importance to accurately predict the bending 
life spans of the electric wires and the wire protecting members . 

Generally, the prediction of the bending life spans of 
the electric wires and the wire protecting members is carried 
out by repeating design, trial manufacturing, and endurance 
test. For example, in such an endurance test, wires distributed 
along an envisioned wiring path are placed on a vibration base 
plate. The vibration base plate is vibrated at a predetermined 
frequency and a predetermined amplitude by the vibrator, and 
a bending life span of the wire when it is vibrated at a specific 
number of vibrations is predicted. 

The literature cited in the present specification are: 
Non-patent document 1 

"Matrix Infinite Element Method" by B. Nasu, published 
by Brain Library Publishing Corp, August 10, 1978, pp 7 to 15 
Non-patent docximent 2 

"Mode Analysis and Dynamic Design" by Kimihiko Yasuda, 
published by Corona Corp, November 10, 1993, pp 54 to 56 

On the other hand, a development period of the vehicle 
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tends to be shorter, and further there is a demand of enhancing 
a bending life prediction accuracy. In this circumstance, the 
related bending life span predicting method in which the 
endurance test is repeated cannot satisfactorily reduce the 
5 development period, and cannot sufficiently meet the demand 
of enhancing the prediction accuracy. 

Disclosure of the Invention 

Accordingly, to cope with the above-mentioned problems, 
10 an object of the present invention is to provide a method and 
apparatus for predicting bending life spans of electric wires, 
which can meet the demands of a highly accurate prediction of 
the bending life span and of reducing the development period 
in an environment where vibrations occur, and a recording medium 
15 storing a program for predicting the same in such ways. 

In order to solve the aforesaid object, the invention 
is characterized by having the following arrangement. 

(1) A bending life predicting method of predicting bending 
life span of a plurality of wires induced by vibration, at least 
20 two points of each of the plurality of wires being constrained, 
the method comprising the steps of: 

a pre-storing step of pre-storing each predicting 
function representing relationships among atmosphere 
temperatures, stresses andbending life spans for the plurality 
25 of wires; 
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a setting step of setting the plurality of wires / the 
atmosphere temperatures/ pre-vibration shapes of the plurality 
of wires / and constraint conditions of the plurality of wires; 
a finite element model forming step of forming finite 
5 element models of the plurality of wires by using a finite element 
method; 

a vibration analyzing step of calculating natural 
frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
10 models which correspond to the natural f recjuencies^ 
respectively; 

a maximum stress retrieving step of retrieving a maximiam 
stress from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires; 
15 a predicting function readout step of reading predicting 

functions corresponding to the atmosphere temperatures set in 
the setting step, respectively; 

a bending life predicting step of acquiring a bending 
life span corresponding to the maximina stress of each of the 
20 plurality of wires while referring to the predicting functions 
read out in the predicting function readout step, and obtaining 
a shortest bending life span from the bending life spans; and 
an output step of outputting the shortest bending life 
span obtained in the bending life predicting step. 

25 

4 
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(2) The method according to (1), wherein in the vibration 
analyzing step, the plurality of wires are regarded as a wiring 
structure in which the plurality of wires are bundled, and 
natural frequencies of the plurality of wires are computed, 

5 respectively. 

(3) The method according to (1), wherein the plurality of 
wires are bundled into a single bundle, and the bundling of 
theplurality of wires is set as one of the constraint conditions 

10 in the setting step. 

(4) The method according to (1) , further comprising a position 
specifying step of specifying a position on the wire 
corresponding to the shortest bending life span, the output 

15 step outputting the position specified by the position 
specifying step. 

(5) The method according to (1), wherein 

in the vibration analyzing step, displacements of finite 
20 elements of the finite element models which correspond to the 
natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
which is induced by vibrations based on the calculated 
25 displacements, and 
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the output step outputs the predicted interference part. 

(6) The method according to ( 1 ) , wherein a curve representing 
a lower confidence interval to a population regression function 
5 statistically calculated using the stresses and data on bending 
endurance life spans that are obtained under a plurality of 
typical atmosphere temperatures for the plurality of wires, 
is used for the predicting function. 

10 (7) A bending life predicting method of predicting bending 
life spans of a plurality of electric wires and wire protecting 
members for protecting the plurality of wires from being bent 
induced by vibrations, at least two points of each of the 
plurality of wires being constrained, the method comprising: 

15 a pre-storing step of pre-storing relationships among 

atmosphere temperatures, stresses and bending life spans for 
plurality of wires and the wire protecting member; 

a setting step of setting the plurality of wires, the 
wire protecting member, the atmosphere temperatures, 

20 pre-vibration shapes of the plurality of wires and the wire 
protecting member, and constraint conditions of the plurality 
of wires and the wire protecting member; 

a finite element model forming step of forming finite 
element models of the plurality of wires and the wire protecting 

25 member by using a finite element method; 
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a vibration analyzing step of calculating natural 
frequencies for the pre-vibration shapes of the plurality of 
wire and the wire protecting member, and calculating stresses 
in individual finite elements of the finite element models which 
5 correspond to the natural frequencies, respectively; 

a maximum stress retrieving step of retrieving maximum 
stresses from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires and the wire protecting 
member; 

10 a predicting function readout step of reading predicting 

functions corresponding to atmosphere temperatures set in the 
setting step; 

a bending life predicting step of accjuiring bending life 
spans corresponding to the maximiam stresses of the plurality 
15 of wires and the wire protecting member, while referring to 
the predicting functions read out in the predicting function 
readout step, respectively, and obtaining a shortest bending 
life span from said bending life spans; and 

an output step of outputting the shortest bending life 
20 span obtained in the bending life predicting step. 

(8) The method according to (7), wherein in the vibration 
analyzing step, the plurality of wires are regarded as a wiring 
structure in which the plurality of wires are bundled, and 
25 natural frequencies of the plurality of wires are computed. 
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respectively - 

(9) The method according to (7), wherein the plurality of 
wires are bundled into a single bundle, and the b\indling of 

5 the plurality of wires is set as one of the constraint conditions 
in the setting step. 

(10) The method according to (7), further comprising a position 
specifying step of specifying a position on the wire or the 

10 wire protecting member corresponding to the shortest bending 
life span, the output step outputting the position specified 
by the position specifying step. 

(11) The method according to (7), wherein 

15 in the vibration analyzing step, displacements of finite 

elements of the finite element models which correspond to the 
natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
20 which or the wire protecting member is induced by vibrations 
based on the calculated displacements, and 

the output step outputs the predicted interference part . 



25 



(12) The method according to (7) , wherein a curve representing 

a lower confidence interval to a population regression function 

8 
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statistically calculated using the stresses and data on bending 
endurance life spans that are obtained under a plurality of 
typical atmosphere temperatures for the plurality of wires and 
the wire protecting member, is used for the predicting function . 

5 

(13) A bending life predicting device for predicting bending 
life spans of a plurality of wires induced by vibrations, at 
least two points of each of the plurality of wires being 
constrained, the device comprising: 
10 a pre-storing unit for pre-storing each predicting 

function representing relationships among atmosphere 
temperatures, stresses and bending life spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 
15 atmosphere temperatures, pre-vibration shapes of the plurality 
of wires, and constraint conditions of the plurality of wires; 

a finite element model forming unit for forming finite 
element models of the plurality of wires by using a finite element 
method; 

20 a vibration analyzing unit for calculating natural 

frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

25 a maximum stress retrieving unit for retrieving a maximina 

9 
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stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for reading predicting 
functions corresponding to the atmosphere temperatures set in 
5 the setting unit, respectively; 

a bending life predicting unit for acquiring a bending 
life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
read out in the predicting function readout unit, and obtaining 
10 a shortest bending life span from the bending life spans; and 
an output unit for outputting the shortest bending life 
span obtained in the bending life predicting unit. 

(14) The device- according to (13), wherein 
15 the vibration analyzing unit calculates displacements 

of finite elements of the finite element models which correspond 
to the natural frequencies, 

the device includes an interference part predicting unit 
for predicting an interference part on the plurality of wires 
20 which or the wire protecting member is induced by vibrations 
based on the calculated displacements, and 

the output unit outputs the predicted interference part. 

(15) A computer readable recording medium storing a program 

25 for predicting bending life spans of a plurality of wires induced 

10 
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by vibrations, at least two points of each of the plurality 
of wires being constrained, the program causing a computer to 
function as: 

a pre-storing unit for pre-storing each predicting 
5 function representing relationships among atmosphere 
temperatures, stresses andbending life spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 
atmosphere temperatures, pre-vibration shapes of the plurality 
10 of wires, and constraint conditions of the plurality of wires; 

a finite element model foiiaing unit for forming finite 
element models of the plurality of wires by using a finite element 
method; 

a vibration analyzing unit for calculating natural 
15 frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

a maximum stress retrieving unit for retrieving a maximum 
20 stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for reading predicting 
functions corresponding to the atmosphere temperatures set in 
the setting unit, respectively; 
25 a bending life predicting unit for acquiring a bending 
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life span corresponding to the maximuun stress of each of the 
plurality of wires while referring to the predicting ftmctions 
read out in the predicting function readout unit, and obtaining 
a shortest bending life span from the bending life spans; and 
5 an output unit for outputting the shortest bending life 

span obtained in the bending life predicting unit. 

Brief Description of the Drawings 

Fig. 1 is a block diagram showing a basic configuration 
10 of the present invention. 

Fig. 2A is a diagram showing an electric wire handled 
in the invention when it is discretized. Fig. 2B is a diagram 
showing beam elements and nodes of the discrete electric wire 
shown in Fig. 2A. 
15 Fig. 3 is a block diagram showing a hardware configuration 

constructed according to the invention. 

Figs. 4A and 4B are diagrams showing bending life span 
data stored in a storing device in Fig. 3. 

Fig. 5 is a flow chart showing a main processing procedure 
20 according to an embodiment of the invention. 

Fig. 6 is a flow chart showing a vibration analysis, 
processing procedure shown in Fig. 5. 

Figs . 7A, 7B and 7C are diagrams showing predictive shapes 
of the electric wire in primary, secondary and tertiary vibration 
25 modes, respectively. 
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Best Mode for Carrying out the Invention 

The preferred embodiment according to the present 
invention will be described with reference to the accompanying 
5 drawings. Pre-conditions, theories and basic formulae which 
are used in the embodiment will first be described with reference 
to Figs. 2A and 2B. Fig. 2A is a diagram showing an electric 
wire in a state it is discretized, handled in the invention. 
Fig. 23 is a diagram showing beam elements and nodes of the 
10 discrete electric wire shown in Fig. 2A. 

In the invention, a finite element method is utilized 
on the basis of the following assumptions: 

(1) The electric wire consists of an elastic member. 

(2) The electric wire is an aggregation of beam elements 
15 serially connected to one another. 

(3) A linearity is retained in each beam element. 

(4) A cross section of the electric wire is circular. 

Next, this electric wire 1 is discretized. Specifically, 
20 as shown in Fig. 2A, the electric wire 1 is divided (discretized) 
into plural beam elements CI, C2, C3, ... . In other words, the 
electric wire 1 is regarded as an aggregation of a finite number 
of beam elements serially connected to one another. 

Accordingly, the electric wire 1 may be considered as 

25 an aggregation of beam elements CI, C2, C3, ... which are serially 

13 
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connected at nodes Nl, N2, N3, .„ to one another^ as shown in 
Fig. 2B. Properties necessary for the beam element are, for 
example, length "1" (see Fig. 2A) , a cross section area A (see 
Fig. 2A) , a moment of inertia of area, a polar moment of inertia 
of area, density, amodulus of longitudinal elasticity, amodulus 
of transverse elasticity, and the like. Those properties may 
be obtained in advance by measurement or calculation. It is 
preferable that the properties are formed as a database, and 
available ad libitum. In the specification, the length and 
the cross sectional area are set as geometric properties, and 
the moment of inertia of area, the polar moment of inertia of 
area, the modulus of longitudinal elasticity, and the modulus 
of transverse elasticity are set as material properties. 

As well known. Hook' s Law holds for a displacement of 
a structure not being vibrating by a static force within a range 
of elasticity of the structure. Likewise, Hook's Law holds 
also for the beam elements CI, C2, C3, as is known. The 
conditions of the displacement continuity and the equilibrium 
of forces are satisfied among the adjacent beam elements CI, 
C2, C3, .... By utilizing those facts, those beam elements CI, 
C2, C3, ...are connected and the resultant discrete body (electric 
wire) is expressed by the following expression (1) . 

[K] {x} = {F} ... (1) 

In the above expression, [K] is a global stiffness matrix which 

14 
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are assembled from stiffness matrices of individual beam 
elements. To state it simply, the contents of the global 
stiffness matrix [K] is like an aggregation of the properties. 
Those values are previously set. {x} is a nodal displacement 
vector (or a displacement vector, simply) , and an array of 
displacement components of all the nodes. The component of 
the nodal displacement vector {x} includes a preset vector, 
such as a point of constraint, and an unknown quantity. {F} 
is a net external force vector, and an array of external force 
components at all the nodes. 

In this way, a finite element model of the discrete wire 
having divided into an optional number of beam elements is formed . 
A technique, which is similar to a technique disclosed in 
Japanese patent application No, 2002-279503, filed on September 
25, 2002 by the J^plicant of the present patent application, 
can be applied to the formation of the finite element model. 
The general matrix finite element method is described in, for 
example, the non-patent documents described above. 

In the invention, to predict a bending life span of the 
electric wire induced by vibrations, it is necessary to obtain 
the natural frequencies, displacements, stress and the like. 
A theory and calculation formulae, which are used for obtaining 
the natural f recjuencies, displacements, stress and the like, 
will be briefly describedbelow, while developing the expression 
(1) . 

15 
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To dynamically handle the wire for predicting its bending 
life span, an inertia force [M] {x"} is added to the load term 
in the expression of the static equilibriiam expressed by the 
expression (1), and damping [C] {x' } is taken into account, 
5 and then we have 



[M] {x"} + [C] {xM + [K] {X} = {F (t)} ... (2) 

where 

[M] : inertia matrix 
10 [C] : damping coefficient matrix 

{x' } : first order differential of the displacement vector {x} 
{x"} : second order differential of the displacement vector 

{X} 



15 In a free vibration state, if external force =0 ({F{t) } 

= 0 in the expression (2) ) , and no damping is present [C] {x' } 
= 0 in the expression (2)) , each point of the wire oscillates 
in a simple harmonic motion at a fixed frequency. When the 
following expression (3) is s\ibstituted into the expression 

20 (2), 

{x} = {Xo} sin G)t ... (3) 

then we have 

[M] ({Xo} sin ot)'' + [K] ({xo} sin ot) = 0 
(-0)2 [M] + [K]) {X} sing cot = 0 
25 ([K] - ©2 [M]) {X} = 0 ... (4) 
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From the expression (4), the natural frequency flo and the 
displacement vector {x} as an eigen-vector can be obtained. 
•Even if the displacement vector {x} is multiplied by a fixed 
number, the expression (4) holds. Hence, the displacement 
vector {x} is a relative displacement. Accordingly, a stress 
obtained from the displacement vector {x} is also a relative 
value. An absolute value obtained from the expression (4) is 

only the natural frequency ©. To calculate actual 

displacements and stresses, it is necessary to input a load 

condition of an actual vibration. 

To dynamically handle the wire for its bending life span 

prediction, if the time-dependant load {F (t) } is input to the 

expression (2) , displacements and stresses each for each unit 

time can be calculated, basically. In case where the equation 

is solved by inputting the time-dependant load {F (t) } to the 

expression (2) , an analysis time to reproduce the vibration 

is long in a phenomenon attendant with periodic vibrations, 

and hence, a time to output the analysis result is extremely 

long. To avoid this, a common practice is to convert the 

expression from a time domain to a frequency domain by Fourier 

transformation. In the calculation of the expression in the 

frequency domain, viz., in the frequency response analysis, 

the frequency-dependant load is input to the expression to 

calculate displacement and stress per each frequency. For the 

17 
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frequency-dependant load, a force applied to points of 

constraint on a wire harness or a wire protecting member as 

an object under prediction of the bending life is set at an 

optional value, the force attending envisioned vibrations of 

the engine, vehicle body, and the like- Thereafter, the 

expression is returned from the frequency domain to the time 

domain by inverse Fourier transformation, whereby an amount 

of calculation is reduced. This technique is well known, and 

hence, no further description of it will be given here. 

By utilizing the theories and calculation formula thus 

far described, one can calculate natural frequency, 

displacements and stresses, which will be described later on. 

The general eigenvalue analysis as described above is described 

also in the non-patent dociament 2, for example. 

In the present invention, a vibration analysis of a 

plurality of electric wires including, for example, grommets 

as wire protecting members is also carried out. To apply the 

finite element method to the wire protecting member, the wire 

protecting member is divided into a plurality of triangular 

finite elements, and stresses at those finite elements are 

obtained. A vibration analysis of the wire protecting member, 

which follows the formation of the finite element model 

consisting of such finite elements, is carried out in a manner 

as in the case of the electric wire. 

Next, a hardware configuration for realizing a process 

18 
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procedure constructed according to the present invention will 

be described. Fig. 3 is a block diagram showing a hardware 

configuration constructed according to the invention. Figs. 

4A and 4B are diagrams showing bending life span data stored 

in a storing device in Fig. 3. 

As shown in Fig. 3, in the present invention, a personal 

computer, for example, is used which has a basic configuration 

made up of a micro-computer 51, an input device 52, a display 

device 53, a printing device 54, a storing device 55, a 

communication interface 56, and a read/write device 57. The 

micro- computer 51 includes a CPU (central processing unit) 51a, 

a ROM 51b for storing a boot program and others, and a RAM 51c 

for temporarily storing various processing results . The input 

device 52 is a keyboard, a computer mouse and the like; the 

display device 53 is an LCD, a CRT and the like for displaying 

processing results; and the printing device 54 is a printer 

for printing processing results. 

The storing device 55 is, for example, a hard disc drive, 

and the communication interface 56 is, for example, a modem 

board for communicating with external devices by way of an 

internet, a LAN line and the like. The read/write device 57 

reads a bending life span predicting program 59a constructed 

according to the invention, which is stored in a recording medium 

59, and writes a result file 55c into the recording medixam 59. 

Those constituent components are interconnected through an 

19 
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internal bus 58 • 

At least a set data file 55a, a bending life data file 
55b, and a result file 55c are stored in the storing device 
55. The set data file 55a contains data preset for predicting 
5 the bending life span. The data represents types of electric 
wires and wire protecting members, atmosphere temperatures, 
pre-vibration shapes, constraint conditions and others, which 
are set in a step SI to be described later. The set data file 
55a contains position information of obstacles 2, such as 
10 pre-estimated stays and electric parts (see Fig. 7) . 

The bending life data file 55b is an aggregation of 
predicting functions yl, y2 and y3 as statistically calculated 
using data of stresses and the number of bendings for endurance 
on electric wires 55al and 55a2, and wire protecting member 
15 (e.g., grommet) 55a3, which are gathered at different atmosphere 
temperatures of, for example, - 40''C, O^C and 25 °C. 

A curve y22 represents an upper confidence interval and 
the curve y23 represents a lower confidence interval to a 
population regression function y21 obtained by the known 
20 regression analysis, as shown in Fig. 4B. For the predicting 
function, the curve y23 of those curves y22 and y23 is used. 
The confidence interval is, for example, 95%. Such predicting 
foonctions are obtained under condition that the electric wire 
and the grommet are placed at different atmosphere temperatures . 
25 Accordingly, thebending life span is predicted under more severe 
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statistical conditions. As a matter of course, a given 
statistical reliability on the predicting function is taken 
into account/ and a computing process of the predicting function 
is easy. As a result, thebending life span is severely predicted 
5 without addition of complicated processing procedures. This 
contributes to further enhancement of quality and wiring path 
plan. Incidentally, the number of bendings for endurance may 
be obtained using the population aggression function. 

The result file 55c contains all the stresses of the 
10 discrete wire and the discrete grommet in the finite elements, 
which are recorded by each natural frequency. The result file 
55c is stored in a text format, and may be output ad libitum. 
The storing device 55 corresponds to predicting function storing 
means set forth in claims. 
15 With such a configuration, the micro- computer 51 installs 

a bending life span predicting program 59a that is read by the 
read/write device 57 into the storing device 55. After power 
on, the micro-computer 51 is activated according to a boot 
program that is stored in the ROM 51b, and starts the bending 
20 life span predicting program 59a. According to the bending 
life spanpredicting program 59a, the micro-computer 51 predicts 
a bending life span of the electric wire and/or the wire 
protecting member induced by vibrations, causes the display 
device 53 and the printing device 54 to display and print the 
25 prediction result, and stores the result in the storing device 
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55. The bending life span predicting program 59a may be 
installed to another personal computer having the same 
configuration, and after installing, operates the computer as 
a bending life predicting device. The bending life span 
5 predicting program 59a may be provided through a cominunication 
line, such as an Internet, a LAN or the like, instead of the 
recording medium 59. 

Processing procedures according to the present 
embodiment of the invention will be described using an 
10 explanatory diagram of Figs. 7A through 7C and flow charts of 
Figs. 5 and 6. Fig. 5 is a flow chart showing a main processing 
procedure according to an embodiment of the invention. Fig. 
6 is a flow chart showing a vibration analysis processing 
procedure shown in Fig. 5. Figs. 7A, 7B and 7C are diagrams 
15 showing predictive shapes of the electric wire in primary, 
secondary and tertiary vibration modes, respectively. 

As shown in Fig. 5, data necessary for the prediction 
is set in a step SI. Specifically, items to be set are at least 
a plurality of electric wires (exactly types of electric wires) 
20 whose bending life spans are to be predicted, wire protecting 
members, (sometimes not recjuired) , atmosphere temperatures, 
pre-vibration shapes for each wires, and constraint conditions 
for each electric wires . To set those necessary items, an input 
screen (not shown) is presented on the display device 53, and 
25 necessary items are entered to the input screen by use of the 
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input device 52. 

The types of the electric wires are the types of the 

electric wires whose bending life spans are to be predicted. 

For the electric wires / geometrical properties and material 

5 properties are linked to the types of the electric wires. 

Specifically, the geometric properties include information 

about lengths and cross sections for each wire. The material 

properties are a moment of inertia of area, a polar moment of 

inertia of area, density a modulus of longitudinal elasticity, 

10 and a modulus of transverse elasticity. Those materials may 

be obtained in advance by tests, for example . Those properties 

concerns the finite elements in the stiffness matrix [K] of 

the expression (1) . The atmosphere temperatures are 

temperatures around the electric wires and/or wire protecting 

15 members when the bending life span prediction is carried out, 

and typically -40°C, 0°C and 25°C. 

The pre-vibration shape is a shape of the electric wire 

or the like in a static state as indicated by reference niameral 

la in Fig. 7A. In the illustration, it is simplifiedand depicted 

20 linearly. The pre-vibration shape la may be manually entered 

by use of the input device 52, or a pre-vibration shape that 

was previously formed before the bending life span prediction 

may be used for the same. Preferably, such a pre-vibration 

shape as to satisfy given constraint conditions and the minimum 

25 bending radius of the electric wire is calculated in advance, 

23 



wo 2004/048940 



PCT/JP2003/015015 



and it is used for the pre-vibration shape la. 

For the constraint conditions, information of perfect 

constraint, rotation constraint and perfect free are set, while 

corresponding to coordinate information on nodes nl to n5 in 

fig 7A. The nodes nl to n5 are connection points of the finite 

elements into which the electric wire is divided. Some of the 

nodes correspond in position to parts at which constraint members, 

such as connectors and clips, respectively. In this instance, 

the nodes nl and n5 are put under the perfect or rotation 

constraint condition by the constraint members, while the 

remaining nodes n2 to n4 are put under the perfect free condition . 

The values set here concern the values of the finite elements 

in the displacement vector {x} in the expression (1) . 

In some cases, the plural electric wires which form a 

wire harness are bundled, by an insulation locking band, tape 

or the like, into a single bundle, and in other cases, those 

electric wires are not bundled. In case where the electric 

wires are bundled together, the wire bundling is set as the 

constraint condition. In setting the constraint condition of 

the wire bundling, it is preferable to handle the electric wires 

that come in contact with the insulation locking band or tape 

in distinction from those located on the inner side of the former 

electric wires. By so doing, the bending life spans of the 

wires can be more exactly predicted since the electric wires 

bundled and the electric wires not bundled are separately 
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processed. 

In some cases, those electric wires are passed through 
a wire protecting member, such as a grommet, and in other cases, 
the electric wires are arranged or distributed without using 
5 the wire protecting member. The plural electric wires passed 
through the wire protecting member are used such that as 
described above, in some instances, those wires are bundled 
together into a single bundle, and in other instances, those 
wires are not bundled • In case where the wire protecting member 

10 is used, a constraint condition suitable for the wire protecting 
member is to be set . The step SI thus far described corresponds 
to "setting step" and "setting unit" in claims. 

Next, in a step S2, finite element models of the electric 
wires and/or wire protecting members are formed. The finite 

15 element model of the wire protecting member is formed by dividing 
the wire protecting member into plural rectangular finite 
elements. The finite element model of the electric wire is 
formed such that the electric wire is divided into plural beam 
elements. The step S2 corresponds to a finite element model 

20 forming step and a finite element model forming means in claims . 

Next, in a step S3, a vibration analyzing process is 
carried out for each finite element. Specifically, in a step 
S301, on the assumption that a wire harness as a bundle of a 
plurality of electric wires as objects vibrates, the expression 

25 (4) applies to the wire harness to calculate a natural frequency 
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of the wire harness. In this case, the geometric properties 
and the material properties as the individual elements of the 
stiffness matrix [K] are calculated on the assimption that those 
properties are the results of composing the geometric properties 
and the material properties of each electric wire and/or wire 
protecting member. For the natural frequency, its values 
corresponding to, for example, the first, second and tertiary 
vibration modes are calculated as shown in Figs. 7A to 7C. The 
number of vibration modes to be calculated and the natural 
frequencies are not limited to the illustrated ones. For 
example, the first, second and tertiary vibration modes are 
set as default values, and may be changed by use of the input 
device 52 . 

Following the natural frequency calculation, steps S302 
and S303 are executed in which the natural frequencies are 
applied to an expression obtained by Fourier transforming the 
expression (2) as described above to thereby calculate 
displacements and stresses of all the finite elements for each 
electric wire and/or wire protecting member. In this case, 
information as to whether or not the electric wires are bundled 
into one bundle and information as to whether or not specific 
electric wire or wires come in contact with the insulation 
locking band or the tape, which are set as the constraint 
conditions in the step SI, are incorporated regarding as the 
information concern all of the inertial matrix [M] , the damping 
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coefficient matrix [C] , and the external force {F (t) } in the 

expression (2) . 

In a step S304, the displacements and the stresses 

calculated every each natural frequency are linked with 
(finite elements) information of the position of each electric 

wire and/or wire protecting member as an ob j ect under prediction 

of the bending life, and are written into the result file. Those 

steps S3, S301 to S303 correspond to a vibration analyzing step . 

and vibration analyzing unit in claims. 

Subsequently, in a step S4, a maximum stress is retrieved 

from those stresses recorded in the result file, for each 

electric wire and/or wire protecting member under prediction 

of the bending life. In a step S5, a predicting function to 

the electric wire and/or wire protecting member under prediction 

of the bending life is read out. Specifically, a predicting 

function of the atmosphere temperature as set in the step SI 

to the electric wire and/or wire protecting member under 

prediction of the bending life is read out from the bending 

life data file 55b stored in the storing device 55. The step 

S4 corresponds to a maximum stress retrieving step and a maximum 

stress retrieving unit in claims, and the step S5 corresponds 

to a predicting function readout step and predicting function 

readout unit in claims. 

In the next step S6, each bending life span corresponding 

to maximum stress of each electric wire and/or each wire 
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protecting member, which is obtained in the step S4, is acquired, 
while referring to the predicting functions read out in the 
step S5. From those bending life spans, a shortest bending 
life span is obtained. The shortest bending life span may be, 
5 for example, the nixmber of bendings for endurance obtained from 
the bending life function in Fig. 4, or may be an endurance 
time calculated based on the number of bendings for endurance 
and the natural frequency corresponding to the former. The 
step S6 corresponds to a bending life predicting step and a 
10 bending life predicting unit in claims. 

In the step S7, a position corresponding to the shortest 
bending life span is specified from its finite element. This 
position specifying is useful in locating a part of the wire 
at which the wire is possibly disconnected. In a step S8, the 
15 position specified in the step S7 and the bending life span 
predicted in the step S6 are output to the display device 53. 
Because the position on the wire protecting member and/or the 
electric wire corresponding to the shortest bending life span 
is output in addition to the bending life span, more accurate 
20 bending life prediction is ensured. The step S7 corresponds 
to a position specifying step and position specifying means 
in claims, and the step S8 corresponds an output step and an 
output unit in claims. 

The following procedure may be taken: in a step S9, an 
25 interference part is predicted, and in a step SIO a predicted 
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inference part is output • Specifically, as described in the 

step S304, the stresses calculated for each natural frequency 

are stored in the result file in a state that those stresses 

are linked with information of the positions (finite elements) 

of the electric wires and/or wire protecting members. 

Therefore, predictive shapes canbe generated for each vibration 

mode from the result file, as indicated by lal in Fig. 7A, la2 

in Fig. 7B, and la3 in Fig. 7C. The position information of 

the obstacle 2, as described above, is stored in the storing 

device 55. Therefore, if the position information is composed 

with the predictive shape for each vibration mode, one can see 

that there is a possibility that the electric wire will interfere 

with the obstacle 2 (see Fig. 7C) . On the contrary, it is 

possible to predict a part of the electric wire where it will 

interfere with the obstacle 2. The step S9 corresponds to an 

interference part predicting step in claims. 

If in a step Sll, it is judged that a command of outputting 

the result file is present, the contents of the result file 

55c is output in the text f oanaat in a step S12 . The output 

contents may be displayed by the display device 53 or printed 

on a sheet by the printing device 54. Further, the contents 

to be output may be designated by the input device 52. If the 

contents outputting is not needed, the secjuence of processing 

steps may be ended (N in step Sll) . 

Thus, the embodiment of the invention succeeds in 

29 



wo 2004/048940 



PCT/JP2003/015015 



providing a method and apparatus for predicting bending life 
spans of electric wires or the like, which can more accurately 
predict bending life spans of the wires and can satisfactorily 
meet the demand of enhancing the prediction accuracy in an 
environment where vibrations occur, and a program for predicting 
the same in such ways. In particular, a disconnection life 
of the electric wire, which is due to engine vibration of the 
automobile, can be predicted for a short time and without 
performing the endurance test. This is a very advantageous 
effect. Further, also statistically, the bending life span 
is predicted under severe conditions. Accordingly, the 
embodiment contributes to further enhancement of quality and 
wiring path plan. 

It will be readily understood that the invention is not 
limited to the embodiment above mentioned. For example, the 
wire protecting member is not limited to the grommet . Further, 
the invention is not limited to the inside of the automobile, 
but may be applied to a place in a factory where vibrations 
occur . 

According to the invention, the plural electric wires, 
the atmosphere temperatures, pre- vibration shapes of the plural 
electric wires , and constraint conditions of the plural electric 
wires are set, and finite element models of plural electric 
wires are formed. Natural frequencies for the pre-vibration 
shapes of the plural electric wires are calculated, and stresses 
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in individual finite elements of the finite element models which 

correspond to the natural frequencies are calculated, and a 

maximum stress is retrieved from the calculated stresses for 

each electric wire . Predicting functions corresponding to the 

plural electric wires and the atmosphere temperatures, which 

are set, are read out, A bending life span corresponding to 

the maximum stress of each electric wire, is acquired while 

referring to the predicting functions read out, and a shortest 

bending life span is obtained from the bending life spans, and 

output. Accordingly, the bending life spans of the plural 

electric wires can be accurately predicted without performing 

the endurance test. As a result, there is provided a bending 

life predicting method which can meet the demands of reducing 

the development period and enhancing the prediction accuracy. 

According to the invention, the plural electric wires 

and wire, protecting members, the atmosphere temperatures, 

pre-vibration shapes of the plural electric wires, and 

constraint conditions of the plural electric wires are set, 

and finite element models of plural electric wires are formed. 

Natural f recjuencies for the pre-vibration shapes of the plural 

electric wires are calculated, and stresses in individual finite 

elements of the finite element models which correspond to the 

natural frequencies are calculated, and a maximiom stress is 

retrieved from the calculated stresses for each electric wire. 

Predicting fiinctions corresponding to the plural electric wires 
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and the atmosphere temperatures, which are set, are read out. 

A bending life span corresponding to the maximum stress of each 

electric wire, is acquired while referring to the predicting 

functions read out, and a shortest bending life span is obtained 

5 from the bending life spans, and output. Accordingly, the 

bending life spans of the plural electric wires inclusive of 

the wire protecting member can be accurately predicted without 

performing the endurance test. As a result, there is provided 

a bending life predicting method which can sufficiently meet 

10 the demands of reducing the development period. 

According to the invention, the plural electric wires 

are regarded as a wire-like structure of a bundle of the plural 

electric wires, and natural frequencies of the plural electric 

wires are computed. This feature results in reduction of 

15 calculation amount. Further reduction of the time taken for 

the bending life span prediction is realized. 

According to the invention, when the plural electric wires 

are bundled into a single bundle, and the bundling of the plural 

electric wires is set as one of the constraint conditions in 

20 the setting step. This ensures more exact prediction of the 

bending life span. Actually, insomecases, the plural electric 

wires are bundled into a single bundle, and in other cases, 

those are not bundled. The bending life span can be predicted, 

while those cases are clearly separated from each other. 

25 According to the invention, a position specifying step 
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for specifying a position on the wire protecting meiober or the 
electric wire corresponding to the shortest bending life span, 
is also specified, together with the bending life span, and 
is output. This feature ensures a more accurate prediction 
5 of the bending life span. 

According to the invention, an interference part induced 
by vibrations on each wire protecting member or each plural 
electric wires is also predicted. Accordingly, with this 
feature, optimum wiring paths can be designed free from wire 
10 disconnection by its contact. 

According to the invention, a curve representative of 
a lower confidence interval to a population regression function 
statistically calculated using data of the stresses and the 
number of bendings for endurance under a plurality of typical 
15 atmosphere temperatures for the wire protecting members and 
the electric wires, is used for the predicting function. 
Accordingly, the bending life span is predicted under severe 
statistical conditions. As a matter of course, a given 
statistical reliability on the predicting function is taken 
20 into account, and a computing process of the predicting function 
is easy. As a result, the bending life span is severely predicted 
without addition of complicated processing procedures. This 
contributes to further enhancement of quality and wiring path 
plan. 
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CLAIMS 

!• A bending life predicting method of predicting bending 
life span of a plurality of wires induced by vibration, at least 
two points of each of the plurality of wires being constrained, 
5 the method comprising the steps of: 

a pre-storing step of pre-storing each predicting 
function representing relationships among atmosphere 
temperatures, stresses andbending life spans for the plurality 
of wires; 

10 a setting step of setting the plurality of wires, the 

atmosphere temperatures, pre-vibration shapes of the plurality 
of wires, and constraint conditions of the plurality of wires; 

a finite element model forming step of forming finite 
element models of the plurality of wires by using a finite element 

15 method; 

a vibration analyzing step of calculating natural 
frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
20 respectively; 

a maximiom stress retrieving step of retrieving a maximum 
stress from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires; 

a predicting function readout step of reading predicting 

25 functions corresponding to the atmosphere temperatures set in 
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the setting step, respectively; 

a bending life predicting step of acquiring a bending 
life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
5 read out in the predicting function readout step, and obtaining 
a shortest bending life span from the bending life spans; and 
an output step of outputting the shortest bending life 
span obtained in the bending life predicting step, 

10 2 . The method according to claim 1, wherein in the vibration 
analyzing step, the plurality of wires are regarded as a wiring 
structure in which the plurality of wires are bundled, and 
natural frequencies of the plurality of wires are computed, 
respectively . 

• 15 

3. The method according to claim 1, wherein the plurality 
of wires are bundled into a single bundle, and the bundling 
of the plurality of wires is set as one of the constraint 
conditions in the setting step. 

20 

4. The method according to claim 1, further comprising a 
position specifying step of specifying a position on the wire 
corresponding to the shortest bending life span, the output 
step outputting the position specified by the position 

25 specifying step. 
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5, The method according to claim 1, wherein 

in the vibration analyzing step, displacements of finite 
elements of the finite element models which correspond to the 
5 natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
which is induced by vibrations based on the calculated 
displacements , and 
10 the output step outputs the predicted interference part. 

6. The method according to claim 1, wherein a curve 
representing a lower confidence interval to a population 
regression function statistically calculated using the 

15 stresses and data on bending endurance life spans that are 
obtained under a plurality of typical atmosphere temperatures 
for the plurality of wires, is used for the predicting function. 

7. A bending life predicting method of predicting bending 
20 life spans of a plurality of electric wires and wire protecting 

members for protecting the plurality of wires from being bent 
induced by vibrations, at least two points of each of the 
plurality of wires being constrained, the method comprising: 
a pre-storing step of pre-storing relationships among 
25 atmosphere temperatures^ stresses and bending life spans for 
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plurality of wires and the wire protecting member; 

a setting step of setting the plurality of wires, the 
wire protecting member, the atmosphere temperatures, 
pre-vibration shapes of the plurality of wires and the wire 
5 protecting member, and constraint conditions of the plurality 
of wires and the wire protecting member; 

a finite element model forming step of forming finite 
element models of the plurality of wires and the wire protecting 
member by using a finite element method; 
10 a vibration analyzing step of calculating natural 

frequencies for the pre-vibration shapes of the plurality of 
wire and the wire protecting member, and calculating stresses 
in individual finite elements of the finite element models which 
correspond to the natural frequencies, respectively; 
15 a maximiim stress retrieving step of retrieving maximum 

stresses from the stresses calculated in the vibration analyzing 
step, for each of the plurality of wires and the wire protecting 
member; 

a predicting function readout step of reading predicting 
20 functions corresponding to atmosphere temperatures set in the 
setting step; 

a bending life predicting step of acquiring bending life 
spans corresponding to the maximum stresses of the plurality 
of wires and the wire protecting member, while referring to 
25 the predicting functions read out in the predicting function 
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readout step, respectively/ and obtaining a shortest bending 
life span from said bending life spans; and 

an output step of outputting the shortest bending life 
span obtained in the bending life predicting step, 

5 

8 . The method according to claim 1 , wherein in the vibration 
analyzing step, the plurality of wires are regarded as a wiring 
structure in which the plurality of wires are bundled, and 
natural frequencies of the plurality of wires are computed, 
10 respectively. 

9. The method according to claim 7, wherein the plurality 
of wires are bundled into a single bundle, and the bundling 
of the plurality of wires is set as one of the constraint 

15 conditions in the setting step. 

10. The method according to claim 7, further comprising a 
position specifying step of specifying a position on the wire 
or the wire protecting member corresponding to the shortest 

20 bending life span, the output step outputting the position 
specified by the position specifying step. 

11- The method according to claim 7, wherein 

in the vibration analyzing step, displacements of finite 

25 elements of the finite element models which correspond to the 
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natural frequencies are calculated, 

the method includes an interference part predicting step 
of predicting an interference part on the plurality of wires 
which or the wire protecting member is induced by vibrations 
5 based on the calculated displacements, and 

the output step outputs the predicted interference part . 

12. The method according to claim 7, wherein a curve 
representing a lower confidence interval to a population 

10 regression function statistically calculated using the 
stresses and data on bending endurance life spans that are 
obtained under a plurality of typical atmosphere temperatures 
for the plurality of wires and the wire protecting member, is 
used for the predicting function - 

15 

13, A bending life predicting device for predicting bending 
life spans of a plurality of wires induced by vibrations, at 
least two points of each of the plurality of wires being 
constrained, the device comprising: 

20 a pre-storing unit for pre-storing each predicting 

function representing relationships among atmosphere 
temperatures/ stresses and bending life spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 

25 atmosphere temperatures, pre-vibration shapes of the plurality 
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of wires, and constraint conditions of the plurality of wires; 

a finite element model forming unit for foirming finite 
element models of thepluralityof wires by using a finite element 
method; 

5 a vibration analyzing unit for calculating natural 

frequencies for the pre-vibration shapes and calculating 
stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

10 a maximum stress retrieving unit for retrieving a maximum 

stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for reading predicting 
functions corresponding to the atmosphere temperatures set in 
15 the setting unit, respectively; 

a bending life predicting unit for acquiring a bending 
life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
read out in the predicting function readout unit, and obtaining 
20 a shortest bending life span from the bending life spans; and 
an output unit for outputting the shortest bending life 
span obtained in the bending life predicting unit. 



14. 

25 



The device according to claim 13, wherein 

the vibration analyzing unit calculates displacements 
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of finite elements of the finite element models which correspond 
to the natural frequencies ^ 

the device includes an interference part predicting unit 
for predicting an interference part on the plurality of wires 
which or the wire protecting member is induced by vibrations 
based on the calculated displacements, and 

the output unit outputs the predicted interference part. 

15. A computer readable recording medium storing a program 
for predicting bending life spans of a plurality of wires induced 
by vibrations, at least two points of each of the plurality 
of wires being constrained, the program causing a computer to 
function as : 

a pre-storing unit for pre-storing each predicting 
function representing relationships among atmosphere 
temperatures, stresses andbending life spans for the plurality 
of wires; 

a setting unit for setting the plurality for wires, the 
atmosphere temperatures, pre-vibration shapes of the plurality 
of wires, and constraint conditions of the plurality of wires; 

a finite element model forming unit for forming finite 
element models of the plurality of wires by using a finite element 
method; 

a vibration analyzing unit for calculating natural 

frequencies for the pre-vibration shapes and calculating 
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stresses in individual finite elements of the finite element 
models which correspond to the natural frequencies, 
respectively; 

a maximiam stress retrieving unit for retrieving a maximum 
stress from the stresses calculated in the vibration analyzing 
unit, for each of the plurality of wires; 

a predicting function readout unit for reading predicting 
functions corresponding to the atmosphere temperatures set in 
the setting unit, respectively; 

a bending life predicting unit for acquiring a bending 
life span corresponding to the maximum stress of each of the 
plurality of wires while referring to the predicting functions 
read out in the predicting function readout \init, and obtaining 
a shortest bending life span from the bending life spans; and 

an output unit for outputting the shortest bending life 
span obtained in the bending life predicting unit. 
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